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ABSTRACT 

Several two-equation models have been proposed and tested against benchmark flow's by 
various researchers. For each study, different numerical methods or codes were used to 
obtain the results which were reported to be an improvement over other models. However, 
these comparisons may be overshadowed by the different numerical schemes used to obtain 
the results. With this in mind, several existing two-equation turbulence models, including 
k — e, k — t, k — u and q — u> models, are implemented into a common flow solver code 
for near wall turbulent flows. Calculations were carried out for low Reynolds number, 
two-dimensional, fully developed channel and boundary layer flows. The quality of each 
model is based on several criterion including robustness and accuracy of predicting the 
turbulent quantities. 

1. INTRODUCTION 

The time averaged Navier- Stokes equations have more unknowns than the number of 
equations. In order to solve this closure problem, it is necessary to model the turbu- 
lent stress tensor, TTfuJ, which appears in the time averaged momentum equation. Many 
semi-empirical models have been proposed, each with its own successes and flaws. The 
two-equation model is one of the more popular approaches. In this model, one equation 
related to the turbulent kinetic energy and one related to the turbulence length scale are 
solved along with the time averaged Navier-Stokes equations. 

This paper summarizes two-equation turbulence models (including recently developed 
models) and compares the robustness and accuracy of different models which have ap- 
peared in the literature. For each model, calculations were carried out for two-dimensional, 
fully developed channel and flat plate boundary layer flows. These flows are appealing for 
model testing because they are simple and self-similar, yet demonstrate important features 
of wall bounded turbulent shear flows. In addition, we can compare the results from these 
calculations with Direct Numerical Simulations (DNS). 

There were four types of two-equation models tested in this study: 

1) k — e 

2) k — u 

3) q — u 

4 ) k — r 
where, 
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k = Turbulent Kinetic Energy, 

q = Vk, 

e = Dissipation Rate, 
u oc e/k = Specific Dissipation Rate, 
t oc k/e — Turbulent Time Scale 


A list of the models which were tested are shown in the table below: 


Ch 

Chien 1 1982 

k — 6 

Sh 

Shih 11 1991 

k — e 

LB 

Lam and Bremhorst 6 1981 

k — e 

NH 

Nagano and Hishida 9 1987 

k — e 

NT 

Nagano and Tagawa 10 1990 

k — e 

LS 

Launder and Sharma 7 1974 

k — e 

r jl 

Jones and Launder 3 1973 

k — e 

MK 

Myong and Kasagi 8 1988 

k — e 

YS 

Yang and Shih 17 1991 

k — e 

WI1 

Wilcox 15 1984 

k — uj 

WI2 

Wilcox 16 1991 

k — uj 

SAA 

Speziale, Abid and Anderson 14 1990 

k — r 

Co 

Coakley 2 1983 

q — uj 


The time averaged momentum and continuity equations axe written as: 


dUi 

dxi “° 

( 1 ) 

D U{ d dUi duiiij 1 dp 

Df dxi dxi dxj p dxi 

(2) 

where the Reynolds stress is modeled as: 


.dUi^dUj. 2,. 

u - u < - ^ a x , + 3 kSii ' 

( 3 ) 


From dimensional analysis, the eddy viscosity is: 

Vt = cu'l' 

where u' and l' are the turbulent velocity scale and turbulent length scale. 


2. THE MODEL EQUATIONS 

In a two-equation model, two turbulent quantities(A: — e, k — t , k — iv or q — us ) are used to 
model the eddy viscosity. The turbulent transport equations of these quantities and the 
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eddy viscosity models are written below. The model constants and other model parameters 
may be found in the appendix. 

2.1 The k-e Model 


vt =(See Table 2) 

DA: d u VT.dk. dUi 

D t + ^ W + " ' ~ e + D 

De _ d u v T ^ de 1 x 1 d*7i # e , ^ 

Di ~a^ + ^7 ) a^ 1 ~ c,/l " C2/2 5) + E 


(4a) 

(5a) 

(6a) 


Wall BC: k = U = 0, € — see table 2 

2.2 The k-r Model 


Jfc 

r = ~ 
e 

V T —Cuftikr 

DA: d . VT.dk dUi k 

D t dxi ^ (7k UlU * dxj t 

Dr _ d .. _ 2 . vt . dr dr 

Di dxi <j T 2 dxi t <j T 2 dxi dxi 

2. VT.dk dr t dUi 

+ a & + ^dTitoTi + (Cl ~ 1 h Ui ' li dr i + {C2h ~ 1} 


(46) 

(56) 


( 66 ) 


Wall BC: k = U = r = 0 

2.3 The k-cu Model 

e 


u 


c,k 


k 

vt =— 
u 


Qk d VT.dk dUi 

TvT =^riy u + — )aT ] ~ 




(T k dxi 

vt a 5a; 




— = — RV + — V 

Dt dxi <j w chc* 

= \(V<j - Vjj) 


£\t T 

- C2U 2 - c 2 c 3 u>(2Si itj n jti )t 


(4c) 

(5c) 

(6c) 


Wall BC: A: = £7 = 0, w. 


6v 
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2.4 The q-cu Model 


u> =y, q = Vk 

k 

(4d) 

(5 d) 

( 6d) 

Wall BC: k = U = 0, gf = 0 




u; 


D ? ^ r/ _L ^ A 

K7=7rr[0 + — ) 


Dt 

Du; 




dx t 


_. 9 UT duj_. 

D t dxi a w dx t ' 


UiUj dUi 
2 q dxj 


C 1 C tl (^- + 


qoj 

Y 

dUj dUi 


y- 


dxj dxi dxj 


cW 


3. MODEL TESTING 

The momentum, continuity, Reynolds stress, eddy viscosity and turbulent transport equa- 
tions (Equations 1-6) are solved simultaneously in a numerical code. The numerical scheme 
is based on GENMIX, a parabolic code developed by Patanker and Spaulding 13 . The tur- 
bulent transport equations and momentum equation are solved by a space marching finite 
difference method obtained by integrating over control volumes. 

Two dimensional channel flow calculations were made at Re r = 180 and Re T = 395. These 
cases were compared with DNS data of Kim et al 4 . Calculations for the two- dimensional 
flat plate boundary layer flow at Re$ = 1410 were compared with DNS data of Spalart 12 . 
Some flat plate boundary layer comparisons were made between experimental data of 
Klebanoff 5 at Re$ = 7700 and solutions of various models. 

Results from channel flow at Re r = 180 and Re r — 395 appear in Figures 1-6 and Figures 
7-12, respectively. Results from flat plate boundary layer flow appear in Figures 13-24. 

An important criterion for two-equation model comparisons is not only how well the model 
predicts mean velocity and shear stress, but also the turbulent kinetic energy and dissipa- 
tion (or specific dissipation) rate. These predictions should provide appropriate turbulent 
velocity and length scales so that the model can be applied to .more complex flows for which 
a simple mixing length model often fails. Some researchers maintain that it is not critical 
that the turbulent kinetic energy and the turbulent length scale are predicted with great 
accuracy. However, one may imagine that a two-equation model making unreasonable tur- 
bulent velocity and length scale predictions would be very questionable when applied to 
more general flows. A model which accurately predicts the shear stress and mean velocity 
does not imply that it has correctly modeled the turbulent kinetic energy and length scale. 
In fact, if the turbulent kinetic energy is incorrect, then the length scale must also be 
incorrect to compensate for the error in the turbulent kinetic energy. For this case, two 
wrongs are making a right. This warrants some caution when computing flows for other 
geometries. 
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The comparisons made in this study are only for rather simple flows. However, We think 
they are important. Because if a model does not correctly predict a simple flow, it can- 
not, in general, be expected to correctly model a more complicated flow. Although the 
comparisons may be highly subjective, it is clear that the JL, LS, WI1 and WI2 models 
underpredict the near wall turbulent kinetic energy compared to the other models. 

The standard k — e model has been proven to provide good results in the high Reynolds 
number range. It is therefore attractive for a near wall k — e turbulence model to approach 
the standard k — e model away from the wall. The LB, LS and YS models are the only 
k — e models in this study which possess this characteristic. 

Because the boundary layer and channel flows are self-similar, the solutions should be 
independent of the initial conditions. However, some of the models (SAA, Co, and LB) have 
difficulty when the initial conditions contain large gradients. The Co Model is particularly 
dependent on the initial conditions. Even slight perturbations to the initial conditions will 
yield noticeably different solutions with this model. 

JL, LS, WI1 and WI2 are the only models which do not contain y + . Damping functions 
which contain j/ + are not desirable because y + is erroneous near flow separations and not 
well defined near complicated geometries. Unfortunately, these are the same models which 
poorly predict the near wall turbulent quantities. 

The Wilcox models (WI1 and WI2) suffer from a numerically awkward boundary condition 
for u at the wall: 


u 


§_u 

C 2 y 2 


as y 


+ 


We cannot define u? at y + = 0. We have tried two ways to approximate u; as y + approaches 
0. First, we chose a large number for uj wa u and, second, we used an asymptotic u wa ii = 
$ . Test cases showed that the solution does not converge as uj wa ii — + oo or — > 0 for 
either model. In addition, both Wilcox models underpredict the turbulent kinetic energy 
peak value for both boundary layer and channel flows. 


4. CONCLUSION 

In our calculations, k — e models such as Ch, NT, Sh and YS were robust and also gave the 
best predictions of overall turbulent quantities. However they all contain an undesirable 
y + in their damping function. 

To explore the capabilities as well as the deficiencies of these models, further testing of 
these models in more complex flows, such as, flows with adverse pressure gradients is 
needed. 
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APPENDIX: Model Parameters and Damping Functions 


Table 1 : k — e Model Parameters 


Model 

n 

D 

E \ 

Ch 

0 

2vk 


Sh 

r - 05 . . Vt /- l • 

l/ u [l-eip(-s+)] cr* ’-d’J 

0 


LB 

0 

0 

0 

NH 

0 



NT 

0 

0 

0 

JL 

0 

-M%£) 2 

2^t(0) 2 

LS 

0 



MK 

0 


■ 

YS 

0 




Table 2 : k — e Model Parameters 


Model 

e 

T t 

V r r 

BCe w 

Ch 

e 

k 

f 

wgnsm 


Sh 

k i ^ i 

€ ~ u ^f- 

k 

e 


i ^HH^S 

LB 

t 

k 

t 

c,U k i 

i/f-4 

aUT 

NH 

€ 

k 

F 

CJA 

0 

NT 

e 

k 

e 

C,U k T 


JL 

e 

k 

£ 

r f * 2 

0 

LS 

e 

k 

£ 


0 

MK 

e 

k 

£ 

C f k * 

U dv 2 

YS 

€ 

7 + (t )’ /2 

CpfukTt 
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Table 3 : k — e Model Parameters 


Model 


Ci 

C2 

V k 


Ch 

.09 

1.35 

1.8 

1.0 

1.3 

Sh 

.09 

1.45 


1.3 

1.3 

LB 

.09 

1.44 

1.92 

1.0 

1.3 

NH 

.09 

1.45 

1.9 

1.0 

1.3 

NT 

.09 

1.45 

1.9 

1.4 1 

1.3 

JL 

.09 

1.45 

2.0 

1.0 

1.3 

LS 

.09 

1.44 

1.92 

1.0 

1.3 

MK 

.09 

1.4 

1.8 

1.4 

1.3 

YS 

.09 

1.44 

1.92 

1.0 

1.3 


Table 4 : k — e Damping Functions 


Model 

fn 



Ch 

1 — exp( — .0115y + ) 


1-.22 eip(--i) 

Sh 

1 — exp(— ,006y + — 4e _4 y + + 
2.5e _e y +3 — 4e~ 9 y +4 ) 

1 

1 - .22 exp(-&) 

LB 

( 1 _ e -.0165«* ) 2 (1 + M5) 

1 + (f) 3 

1 - exp(-Rf) 

NH 

[1 - exp{-^)) 2 

1 

1 — .3 exp(—R 2 ) 

NT 

[l-exp(-£)] 2 (l + ^) 

1 

(1-.3 e*p(— (^) 2 )) 
(1 - exp(=j£)) 2 

JL 

eX P( l + Rjso) 

1 

1 — .3 exp(—R 2 ) 

LS 

eX P((l+fl,/50p) 

1 

1 — .3 exp(-Rj) 

MK 

(* + ^X 1 “ eX P(~ 70 )) 

1 

(1-f exp(-f)) 
(1 -exp(^)) 2 

YS 

1 — exp{— .004y + — 5e _5 y + + 
2 e _6 y+ 3 - 8e _8 y +4 ) 

1 

l-.22exp(-§) 


Table 5 : k — u> Model Parameters 


Model 

c. 

Ca 

c 2 

c 3 

<7* 

C 

WI1 

.09 

5 

. ... ,9 

40 

0 

2.0 

2.0 

WI2 

.09 

<tic 2 (l + v/C„) - Jj-cVCJ 

3 

40 

1 

2.0 

2.0 
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Table 6 : k — t Model Parameters 


Model 

°k 

<T r j 

cr r2 

Ci 

C 2 

c, 

SAA 

1.36 

1.36 

1.36 

1.44 

1.83[1 jfexpt-sjJ-)] 

.09 


Table 7 : k — t Model Damping Functions 


Model 

h 

h 

SAA 

l-expi^) 

(- 1 - + ~ ex P(^70 ')] 


Table 8 : q — u> Model Parameters 


Model 


<J\jj 

c. 

Ci 

C 2 

/, , 

Co 

1.0 

1.3 

.09 

.405/„ + .045 

.92 

1 — exp(— .00657?*) 


R k = Vkl , 
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FIGURE 5 

2-D Channel Flow Re^. = 180 
a ) ^wall — 10 3 , b) W wa ii = 10 5 , c) inwall = 10 7 , d) 


























FIGURE 11 _ 

2-D Channel Flow Re r = 395 
a ) W wa ll — 10 3 , b) OJ wa ll — IQ 5 , c) Uf wa u = 10 7 
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FIGURE 18 

2-D Boundary Layer Flow Reg = 1410 
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FIGURE 20 

2-D Boundary Layer Flow 
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FIGURE 22 
2-D Boundary Layer Flow 

















REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-01QQ 


Public reporting burden for this i ooltecbon of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources 
garnering and maintaining the data needed, and completing and reviewing the collection of Information, Send comments regarding this burden estimate or any other aspect of this 
a>ltectton of Information Jr^uding suggests for reducing this burden, to Washington Headquarters Services, Directorate tor information Operations and Reports 1 21 ^Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704*0188), Washington^ DC 20503 


1. AGENCY USE ONLY ( Leave blank) 


TITLE AND SUBTITLE 


2. REPORT DATE 

September 1991 


A Critical Comparison of Two-Equation Turbulence Models 


REPORT TYPE AND DATES COVERED 

Technical Memorandum 

~ _ I 5 FUNDiNG LUMBERS 


6. AUTHOR(S) 

N.J. Lang and T.H, Shih 


WU-505-62-21 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 

E-6562 


9. SPONSORING/MONITORING AGENCY NAMES(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


10, SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA TM- 105237 

ICOMP-91-15 

CMOTT-91-05 


11. SUPPLEMENTARY NOTES ~~ 

N.J. Lang and T.H. Shih, Institute for Computational Mechanics in Propulsion and Center for Modeling of Turbulence and 
Transition, NASA Lewis Research Center (work funded under Space Act Agreement C-99066-G). Space Act Monitor 
Louis A. Povinelli, (216) 433-5818. 

12.. DISTRIBUTION/AVAILABILITY STATEMENT " 1 12b. DISTRIBUTION CODE 

Unclassified - Unlimited 
Subject Category 34 


13. ABSTRACT (Maximum 200 words) ~ 

Several two-equation models have been proposed and tested against benchmark flows by various researchers. For each 
study , different numerical methods or codes were used to obtain the results which were reported to be an improvement 
over other models. However, these comparisons may be overshadowed by the different numerical schemes used to obtain 
the results. With this in mind, several existing two-equation turbulence models, including k -e, k -t, k ~to and q -co 
models, are implemented into a common flow solver code for near wall turbulent flows. Calculations were carried out for 
low Reynolds number, two-dimensional, fully developed channel and boundary layer flows. The quality of each model is 
based on several criterion including robustness and accuracy of predicting the turbulent quantities. 


14. SUBJECT TERMS 

Turbulence models 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 
NSN 7540-01-280-5500 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


15. NUMBER OF PAGES 

34 

16. PRICE CODE 

A03 


20. LIMITATION OF ABSTRACT 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 




